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Abstract 

The Wigner function formalism has been applied to the analysis of elastic scatter- 
ing processes. The new element of known formalism is the choice of the phase space 
on which the Wigner function is defined. This phase space is 4-dimensional space of 
canonically conjugate variables - the transverse momentum of scattered particle and 
the projection of the vector of the closest approach of particles in the plane which 
is perpendicular to the momentum of initial particle. The exact expression for the 
mean value of the scattered particles production area radius through a scattering 
amplitude has been obtained. 
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1 Introduction 



Wigner function, introduced in 1932 [T] , defines quasiprobability distribution in phase 
space of canonically conjugate (coordinate-momentum) variables. The main object to 
construct the Wigner function is the wave function of the system. With the help of the 
Wigner function one can compute an average value of any dynamical observables. 

The development of the formalism of the Wigner function considerably widened the 
field of its application. This formalism found an application in such fields as hydrody- 
namics j2], problem of N bodies in Boltzmann-Vlasov equation [3], plasma physics [I]. 
It was applied also to calculations of quantum corrections to the transition coefficients 
[5j, as well as, to quantum optics [SI El E] arid statistical mechanics [9]. In Ref. |10| the 
formalism of the Wigner function was applied to the quark and gluon distributions in the 
proton. In this work the relation of Wigner function to generalized proton distribution 
was demonstrated. 

In the present work we generalize the formalism of the Wigner functions to the case of 
an elastic scattering of particles. We suggest to use, as a building block for a construction 
of the Wigner function, a scattering amplitude instead of a wave function of the system. 
In the case of the scattering amplitudes the connection between the momentum and 
coordinate space is not given by the Fourier transformation, as in the case of the standard 
definition of the Wigner function. The point is that the momrntum space in the case 
of scattering processes is not flat, the transformation from the momentum space to the 
canonically conjugate one is performed by the Shapiro function [11]. We shall demonstrate 
that the canonically conjugate space space in the case of the elastic scattering corresponds 
to the space of the nearest approach parameter. Analogous formalism has been developed 
in Ref. |12| in the contex of construction of the relativistic configuration space, as well as 
for the case of Wigner function in the Lobachevsky space |13j . 

The problem of construction of the canonically cojugate space for the case of the elasic 
scattering was solved in Ref. |14] . It was demosnstrated that the transformation from 
the momentum to the "coordinate" space is performed by the Shapiro functions which 
form the complete and orthonormal basis in the two dimensional space of the transverse 
momentum. The Shapiro functions are plane waves on the SO (2.1) group. 

2 The 50^(2.1) algebra formalism 

First we reproduce briefly the formalism of interaction area spatial structure description 
in the framework of the 5*0^(2.1) algebra. The formalism has been described in details 
in the papers pHf [T5| ITSj IT7] . We consider here only elastic two-particle collisions. 

A classical trajectory of free particle escaping from the interaction region continuated 
to the interaction region is located at certain time at a minimum distance from the O 
point (the center of the target, Fig. 1). The coordinates of this point are described by 
the components of the vector d. This vector is expressed through the particle momentum 
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q and its orbital momentum L in the following way: 

di = ~^ £ ijkqjL k , k = 1, 2, 3 . (1) 
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Figure 1: Classical trajectory X(t) of the asymptotically free particle with the momentum q 
continued to the reaction region is characterized by the minimum distance d from the chosen 
point O. This point corresponds to center of the target in the laboratory reference frame and to 
the beam collision point in the center-of-mass system. 



When we consider two-particle elastic processes the vector d represents the vector of 
the closest approach. In the center-of-mass system it defines the mean minimum distance 
between scattered particles. In this case the qi in the ratio (pQ) is relative momentum in 
the center-of-mass system and the Lj is relative orbital momentum. 

The next step is the quantization of d iy it is based on the substitution of C-numbers 
by the corresponding operators in the expression ([T]). The ordering of the operators in 
the expression ([T]) is fixed by the requirement that the operator dj is Hermitian: 

di = ^(EijkqjLk - iqi) , di = {di) + . (2) 

(here and hereinafter H — 1). 

Six generators d±, g^, ds, Li, L2, L3 form a SO (3.1) algebra with the trivial Casimir 
operator. Operators d\, o^, L 3 constitute a nontrivial algebra 50(2.1) with two Casimir 
operators q 2 and K: 

[did 2 ] = ~-^L 3 , [diLz] = -id 2 , [d 2 L 3 ] = %d x . (3) 

K = dl- ^L 2 3 , dl = d 2 + d 2 2 , [L 3 K] ^ , 

[d 1>2 K] = , [d 1(2 q 2 ] = . 

Wave function of a state with fixed value of the Casimir operator and squared mo- 
mentum satisfies the equations 

^-l^W&V, q 2 iP = q 2 iP. (5) 
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We interpret the eigenvalue b 2 , which belongs continuous spectrum, as an escape pa- 
rameter for a particle scattering process on the external field and as a squared vector of 
the closest approach in the case of two-particle elastic scattering. 

Eq. (0) in the momentum representation on the surface q 2 = const are realized with 
the help of the following operators: 



l ( , 2 2\ & o \ 

d2 = 7^ {q2 ~ ql) W.~ qiq2 -k~ q ^ 



t .03 
\ dq 2 oqi 

Solutions of the system flSJ) are the plane waves on the SO (2.1) group [TT| [12} [T8] : 

ip(q±) = €(q±,P), 
where the plane wave £(q±, p) is given by: 



q q~n-q± 



vV -q±\ vV -q\) (7) 

-» -» 2 7 2 2 ^ 

fi = nfi, n = b q - - , 

here the n is arbitrary unit vector. The set of functions £(<fj_, fx) forms a complete ortonor- 
mal system. These functions allow us to perform transformation from the space of the 
functions defined in the momentum space q± to the space of the functions defined in the 
//-space, where 

q± — (q±cosip, q±sinip), < q±_ < q, < (p < 2n, 

jl = (ll cosip, 11 sinip), 0</i<oo, < if) < 2n . 

Direct and inverse transformations between the function f(q±) in the momentum space 
and its image in the /i-space g(p) are given by: 



f (t) (q±) = J dn^(q±^) g {t \fi , 
g {e) (t) = (2^2 / dJ Ve(g±,$ f i€) (qi 



(8) 



e = ±1 , dfig- = — ^=^= dq± , dQ^ = th^Li) djl 



where 

e = ±1 , dfi ? = 

q\/q" - q± 

The signature e correspondes to the sign of the third component of the momentum 
q^. The transformations generated by the algebra do not change the sign of third 
component of the momentum q 3 . Therefore, division of the momentum space into #3 > 
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and (fe < subspaces is invariant under transformations generated by the d\ , d 2 and L 3 
operators. We introduce the f^' amplitudes with fixed signature: 



f (±) (q±) = f(q±,qs = ±\Jq 2 -q ± 2 ) ■ (9) 
We note that for arbitrary function F(q) the following relation is satisfied: 

J F{q) dq = J q 2 dq dfi F(q) = ^ J q 2 dq dfi 9 - F {e \q ± ) . (10) 



e=±l ' 



Here 



dQ = sinO d9 dip , dQq 



dql . 



q\/q 2 - q 2 ± 

Further we introduce the expression for an elactic cross-section of a two-particle process 
in terms of elements of the S'-matrix. As it is showed in Ref. [2], elastic cross-section has 
the form 



dQ 



k ■ 



'ID 



Here the p is the momentum of the incoming particle A in the center-of-mass system, k 
is the momentum of scattered particle A, and |pj = \k\ = p. 

The amplitude f^- ± \k±_',p) is related to the S'-matrix as follows: 



f {± \k r ,p) = 2np A(p) (k ± , ±y/k 2 - k\- h = -k\A\ p- -p) 

(f\F\in) = 5\P m - P f )(f\A\in), S = I + iF . 
The factor X(p) is related to the flux of colliding particles \u\: 



Hp) = tpi 



Ea ■ E 



B 



p(E A + E B ) 



Detailed derivation of the relation fill I) has been given in Refs. |14[ 



3 The formalism of the 50^(2.1) group in the small 
transverse momentum approximation 

We define the profile function g^(fl) of two-particle state on the 5*0^ (2.1) group for 
elastic processes as amplitude transform from the momentum space to the canonically 
conjugate /2-space, which is given by Eq. (JHJ). This function plays a main role in the 
construction of the Wigner function. We will show that the formalism is simplified in 
the limit of small transverse momentum. This limit is justified in such processes where 
scattering occurs on small angles {9 ~ n 9 w 7r) into the forward or backward hemi- 
spheres. Both these cases correspond to the condition qj_/q 1. It is useful to intriduce 
the momentum variables on the unit hyperboloid: 

2 2 - 2 1 

U = Up. — It = 1 . 
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We have the folowing relations for these variables: 
du dq± du q± u 







vi + 



9 % 



5( 2 )(«-iT) = g 2 (l- (/2 



where v = (v , v) 




(13) 



The equation ([5]) for the eigenvalues of the Casimir operator K in the terms of these 
variables has the following form: 



du\ 



rV> + 6 2 gV = 



Its solution (plane wave on the SOJ2.1) group) is given by: 



ip(u) — Uo (n ■ u) 2 



n 



K l, n) , n 



, n ■ u = UqUq — nu = Uq — nu 



b 2 q 2 



1 



(14) 

(15) 
(16) 



(In the terms of momenta q this solution is presented in Eq. (I7j)). 

Now we consider the solution of the equation (114)) in the small transverse momenta 
approximation, that simplifies expression for ip(u). In turn it leads to essential simplifica- 
tion of expression for the Wigner function and analytic expressions for the profile function 

We look for a solution of Eq. (1 14)) in the following form: 



ip(u) 



(p{u) 



(17) 



Substituting this expression to the equation (jl4j) . we obtain for the <p(u) following differ- 
ential equation: 



K + 1) 



+2u 1 u 2 



dp) d 2 p> 



p + 2z/i x - h 



du\ du 2 



+ {ul + I) 



p> + 2ifi 2 



dp d 2 p 
duo du 



2 J 



+ 



dip . dp d 2 p 

-/il/i 2 V 9 + I'P'l-R 1- *A*2— 1- 



du-. 



' du\ du\du2 



p + b 2 q 2 p> = . 



u 



Our aim is to find <p(u) in the region of \u\ <C 1. In accordance with this we look for 
a solution for <p{u) in the form of series over powers of the vector u: 



(p(u) = 1 + BijUiUj + 0{u A ) 



(19) 
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where the Bij is symmetric matrix. Substituting this expression to the equation (IT%|) . we 
obtain the system: 

Li 2 - b 2 q 2 + 2 - 2Sp B = , 

frBu + fi 2 B 12 = , (20) 
ti\B 12 + ^B 22 = . 



These expressions are obtained by equating coefficients in the (I18j) in front of given power 
of the vector u. 

In what follows we shall require that the matrix B = 0, i.e. the function <p(u) in 
Eq. (TTTj) is unity with the accuracy of 0(u A ). With that accuracy 

= e ipa . (21) 

That implies that the transformation from the momentum space u to the canonically 
conjugate /I-space has simple form of the exponential function. The escape parameter 
squared b 2 in this case is the eigenvalue of the Casimir operator K with the accuracy of 
0(u 2 ): 

Ke 1 ^ = b 2 + O P J ■ (22) 

The following conditions are satisfied: 

b 2 q 2 -2 = fi 2 , Sp B = . (23) 

Note that at this point we redefine the parameter fi as compared to Eq. (I16p . We make 
such redefinition in order to interpret the escape parameter b as the distance between 
scattered particles. 

In this approximation we can interpret amplitude fW(u) representation in the form 

F^\u) = J e^V ±} (/2) d/2 (24) 

as expansion of the amplitude in eigenstates of Casimir operators fl3J), (JSJ) with corre- 
sponding interpretation of the vector ft. Inverse transform 

9 {± \t) = J^yJ e^F^(u) du (25) 

defines the function g^(jT). We call it as the profile function on the 5*0^ (2.1) group in the 
small transverse momenta approximation. The momentum space area u = u(cosip, sirup), 
< u < oo, < ip < 2tt, and the parameter space area fl = fi(cosip, simp), < \i < 
oo, < if} < 2ir, are canonically conjugate to each other. 

Approximate Eq. (1241) . in comparison to the exact equation (jSJ), essentially simplifies 
calculations and retains main features of the quantum character of the parameter b: sig- 
nature e, correct description of small phase spaces bq, i.e. region ju w and its physical 
interpretation as the vector of the closest approach. 
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Normalization of the amplitude F^'(u) is chosen in such way that an elastic cross- 
section is given by: 

[ \F^(u)\ 2 du. (26) 



Then the amplitude F^(u) is connected to the introduced above amplitude f^(q±',p) 
by the following expression: 

Fl±,(S) = WTirw /(± ' = p 7rff? ; pL = n = p ) ■ (27) 

Argument of the function on the right hand side of the (127|) is practically determinated 
by the substitution \cos9\ = —. 

Expressions (124")) allow to construct Wigner function in standard way. The amplitude 
and the profile function g^' play the role of the wave functions of the state in the 
momentum and coordinate spaces accordingly. 



4 Wigner function on the 50^(2.1) group 

Let us introduce the S'0 At (2.1)-transform of the correlator of scattering amplitudes: 
W^(u,fl) = j^y 2 j e^> (±) U- (u + ^j dv , (28) 

where the F^ is amplitude determined by Eq. (I2TI) . This correlator is the generalisation 
of the Wigner function for the case of scattering amplitudes. In Eq. (128)) in contrast 
to usual definition of the Wigner function we use scattering amplitudes instead of wave 
functions. Introduced in that way function (u, jl) has following properties: 

1. W {±) = W {±) - (reality) , 

2. / W (±) (u,jl) dfl= \F {± \u 



I 2 _ d(T el 



du 

(±) (29) 



f 2 da ( ' 

3. / W (± \u,jl) du = (2tt) 2 \g (±) (jl)\ = —fir- , 9 {±) is determined by Eq 
J dfj, 

4. J W {± \u,jl) dpi du = a 



'el 



Thus, this function satisfies all properties of the Wigner function and it allows one to cal- 
culate an average value of any function h(u, jl) ordered according to the Weyl prescription 

(h(u, jl)) {±) = -py / h(u, p) W {± \u,jl) dp du . (30) 

a el 

Let us calculate the (b 2 ) - average value of squared minimum distance between scat- 
tered particles in elastic collisions. In accordance with the (12"3"|) this distance is directly 
connected with the (/x 2 ): 

( M 2 ) = g 2 <6 2 ) - 2 . (31) 



It is follows from Eq. Q3QJ) that 

<^ 2 ) (±) = _L { W {± \u,pi) d/x du . (32) 

Substituting here representation (1281) and taking into account that 



we obtain 



<^} (±) = 4) / |V^ {±) (^)| 2 dn 



(33) 



In the case when the amplitude (u) does not depend on the direction of vector u, one 
has: 

a = «, <^ 

Then representation (1331) is reduced to the following one-dimensional integral: 

du . (35) 



oo 



<9w 



Analogous computations give 



fc*) (±) = -^) / F(±) («) fa F^\u) du , fi^ip fc = 1, 2. (36) 
We note that the expression for (// 2 ) can be put in the form similar to Eq. (136]) : 



y) (±) = ^ I F^(u)^F^(u)du. (37) 



This expression corresponds to standard quantum-mechanical representations. Thus, the 
Wigner function allows to obtain dynamic variables depending on the /2 and u in the 
operator form and in any representation. 

We emphasize that, althought the definition of the Wigner function (|2"8"|) is an ap- 
proximate^, Eqs. (I33|35p are exact in the frameworks of Wigner function formalism. The 
reason is that the scattering amplitude drops fastly if q± — ¥ q (angle 9 f» tt/2). Therefore 
the physical interpretation of the impact parameter /i retains even in the approximation 
we use. 

Calculating (b 2 ) it is important to reveal physical character of the parameter, i.e. its 
connection with dynamic parameters of a system. Let us consider simplest models for 
amplitudes - the one-particle exchange of a particle with the mass m in the t-channel and 
the Regge amplitude with the Pomeron exchange. 



lr The approximation is related to the choice of the approximate eigenfunction of the Casimir operator 

(EH) 
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The amplitude of one-particle exchange in the t-channel is: 

/(?. P) = ■ ~ 2 = 7^ — 2 ' 1^1 = Isl = P • ( 38 ) 

t — m z (p — q) z + m z 

Here the g is the coupling constant, p is the momentum of initial particle, q is the mo- 
mentum of scattered particle. Introducing signature and conical variables uq h u, we 
obtain: 

g/2p 2 



/ (e) (g±; Pi = , p 3 = P) = — , e = ±1 • (39) 

From here we obtain 

-3/2 

F «(iZ) = __ ^ _^ _ . (40) 

zp 1 ~~ ^ + W 
Using this expression, we obtain for the elastic cross-section 

(e) 9 2 27T 

*<* = V z (*o - c) ' (41) 



where 



2 

m 



(If m = rriA + "^c parameter zq corresponds to the large half axis of Leman ellipse). 

Substitution of the amplitude (1401) to the right hand side of Eq. and taking into 
account Eq. (|4l|) leads to the following result: 

(/^ 2 ) (+) = ~ *bfo-l) (~| + 3{zo A _ 1)2 + 2^o + 4zg + (4zg - 2^ ) In ( £ ^)) ■ (43) 

In the region of large energies of colliding particles (m 2 /p 2 <C 1) the second term in the 
brackets of Eq. (H3"j) gives the main contribution. In this approximation we have: 

2 ft 



<" 2 > ,+, = V§-^- (44) 

In the case of scattering into backward hemisphere similar computations lead to following 
result: 

</^ 2 > H = ~ *b(*b + l) (-~ - 2z + Azl + 3{zq \ 1)2 + (2^o - 4^) In (^)) • (45) 
In the large energies region we obtain: 

(b 2 > H = 0.39 % . (46) 

Thus, from Eq. (1441) one can see that in the case of one-particle exchange average "im- 
pact" radius of scattered particles into forward hemisphere is determined by the Compton 
length of exchange particle (virtual "cloud"). It is a well known fact and its confirmation 
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allows us to use the expression ( |4*4|) for calculations with more complicated amplitude 
models F(u). Elastic scattering into backward hemisphere (reflex scattering) is deter- 
mined by the impact radius which drops with the increase of energy (146 p . The existence 
of non-trivial backward scattering is a consequence of quantum character of the closest 
approach vector, and it is an essential distinction from eikonal models. As it has been 
shown in Eq. [IT], correct using unitarity condition leads to an important contribution of 
scattering into backward hemisphere. 

As a second example let us consider the model of the elastic pp-scattering correspond- 
ing to exchange of the Pomeranchuk vacuum pole in the t-channel. In this case the 
amplitude of the elastic scattering is: 



f(cos9) = t(/(.s)< 
B = a'(0) 
t=-2p\ 



Bt 



.7T S 

i h In — 

2 s 



(47) 



cost 



The phenomenological parameter x determines slope of the diffractive cone in differ- 
ential cross-section of the pp-scattering, s as is invariant mass, So ~ 100 GeV 2 determines 
the border of the high-energy regime area. The parameter a (0) is the slope of the Pomer- 
anchuk linear trajectory: a(t) = 1 + at (0)t. According to ISR and SPS experimental data 
at energies y/s « 540 GeV the parameter x is x = 12-M6, a (0) = 0.12±0.03 GeV~ 2 [21]. 

Transformation to cone variables gives for the amplitude F^'Ou): 



Taking into account Eq. 



2p 2 B^- 



-3/2 



-2p 2 B 



X = ig(s)e 

and the relation du = u du d(p, we obtain: 



(48 



a 



(p) = 2tt|A| 



du 



II a 



o a/u 



2tt| A| 



(49) 



where a = 4p 2 e(Re B). Elastic cross-section of scattering into forward hemisphere is: 



a 



(+) 

el 



(P) 



Ti g 2 (s) 



2p 2 a'(0) 

Taking into account Eq. fl35j) . we obtain for (fi 2 )^ 



In — + x 

so 



(50) 



2\(e) 



2tt 



IAI 



Uq 



d 



Our 



U 



-3/2 



2p z tB 
' "0 



(51) 



In the high energy approximation we have \B\ 1 and a ^> 1. In this approximation we 
obtain for the forward hemisphere: 



2\(+) 



V 2 ) 



7r 2 /4+ fx + ln — 
2pV(0)- V 



x + In ■ 



(52) 
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In accordance with ( l3Tj) the average value of the squared minimum distance between 
scattered particles is given by: 

7r 2 /4 + (x + ln^ N 



(& 2 ) {+) = 2a'(0) 



x + In 



(53) 



This distance is determined by the parameter of diffractive cone slope x, the Pomeranchuk 
trajectory slope a (0) and by the energy factor lns/so corresponding to Regge-exchange 
dynamic. 

At s ~ So the main contribution to (b 2 )^ is determined by the parameter x. For this 
case we have 

(0.34 



'(b 2 } {+) 

At the SPS-collider energy yfs « 540 GeV: 



(b 2 ) {+) 

At the LHC energy yfi « 14 Tel/: 



(0.44 



0.39) /m . 
0.48) /m . 
0.54) fm . 



(54) 
(55) 
(56) 



'(b 2 ) {+) « (0.50 

For reliable estimations of the (b 2 )^ it is necessary to use realistic amplitude F(u) of the 
elastic process. The amplitude should be in agreement with contemporary experimental 
data and it should satisfy unitarity condition. This condition controls contribution of 
the many-particle intermediate states to the elastic process and it gives considerable 
contribution to (b 2 )^ as energy function. 

Further we calculate the Wigner function W^ + \u, fx) for the model of one-particle 
exchange ( 138]) . For its calculation we should note that the amplitude of scattering into 
forward hemisphere (l4"0"j) is simplified in the following way: 

; 3/2 _ -g/p 2 



F (+)( 



9 Up 
2p 2 + 



u 2 + ^ 

p 2 



(57) 



In this case we obtain for the Wigner function: 



W i+ \u, fx) 



9 Z 



(2tt) 2 p 4 



dA 



u + 



-r p 2 



u 



+ 



(5J 



Reality of the W^ + \u, fx) is ensured by the symmetry of the denominator under transfor- 
mation A — > — A. For the calculation of the integral ( 15 8 j) we use the following expression: 



da 



I)' + (!-«) f« + |) 2 + $ 



(59) 
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Then we have: 

W i+) (u, p) 

9 2 



(2tt) 2 p 4 



Re / dA e i>1A 



da 



a(u-%) + {l-a)(u+§) 



2 • 



(60) 



The operation Re ensures the reality of the function W^ + \u, ft). Further we make the 
following substitution: 



5= - + (l-2a)u 



and obtain: 

W {+) (u, jt) 

4g 2 



(2tt) 2 p 4 
Using the expression 



Re da / d5 e -*Wl-*x) _ 



5 2 + f + 4au 2 (l-a) 



dt t e-^ 2+ ^ 



(5 2 + a 2 

u 

where we introduce a 2 = m 2 /p 2 + 4au 2 (l — a), we obtain: 

1 oo 

±g 2 



W {+ \u, jt) 



(2vr) 2 p 4 



Re / da e - 2 ^(i- 2 °0 / dt t e~ ta / d5 e^ - * 



1 oo 

(2tt) 2 p 4 J J 



ta 2 -n 2 /t 



The integral over variable t is expressed through the cylindrical function |22| : 

v/2 



dx ■ x v-l e -P*~l/* = 2 



K v {7.yff¥f 



Then 



W (+) (u, ft) 



±g 2 



2vr p 4 



da 



\jt\ cos [2/H(l — 2a)] 
Aa(l-a)u 2 + ^- 



K x 2|/xU 4a(l -a)u 2 + — 



(61) 



(62) 



(63) 



(64) 



(65) 



(66) 



This expression is explicit form of the Wigner function for the model of the one-particle 
exchange. This function of the \ft\ h \u\ at different angles ip between these vectors is 
presented in Figs. 2 and 3 (with arbitrary normalisation). 

One can see from Eq. f )66|) and numerical calculation that the Wigner function is 
sign-changing function. It is its usual property. 
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Figure 2: The Wigner function for the model of the one-particle exchange. Dependence on |/2j 
and \u\ at the ip = is shown. 



Figure 3: The Wigner function for the model of the one-particle exchange. Dependence on |/2| 
and \u\ at the ip = tt/4 is shown. 

5 Conclusion 

In the present work we constructed the generalization of the Wigner function, using 
as the building blocks the elastic scattering amplitudes. The canonically conjugate (to 
the transverse momentum) space is formed by the vector of the nearest approach of the 
scattered particles. Constructed in this way the Wigner function allows one to compute 
space-time characteristics of the scattering processes such as the mean radius of the par- 
ticle production region, multiplicities of produced particles with fixed production radius. 
The developed here approach also allows one to obtain the tomographic image of the 
particle production region with help of Radon transformation, this formalism we shall 
discuss elsewhere. 
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